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Abstract. We study deeply virtual Compton scattering and deep exclusive meson electroproduction on 
a deuteron target. We model the Generalized Quark Distributions in the deuteron by using the impulse 
approximation for the lowest Fock-space state on the light-cone. We study the properties of the resulting 
GPDs, and verify that sum rules violations are quite small in the impulse approximation. Numerical 
predictions are given for the unpolarized cross sections and polarization asymmetries for the kinematical 
regimes relevant for JLab experiments and for HERMES at HERA. We conclude that the signal of coherent 
scattering on the deuteron is comparable to the one on the proton at least for low momentum transfer, 
providing support to the feasibility of the experiments. The short distance structure of the deuteron may 
thus be scrutinized in the near future. 

PACS. 2 4.85. +p, 12.38. Bx, 25.30.-e 



1 Introduction 

The study of hard exclusive processes, such as deeply Vir- 
tual Compton Scattering (DVCS) 

eA -> e'jA' (1) 

and deep exclusive meson electroproduction (DEMP) 

eA -> e'MA' (2) 

where A is a hadron (usually a nucleon, here a deuteron), 
and M a meson (usually a p or a 7r) or a pair of mesons (of 
relatively small invariant mass) in the kinematical domain 
of a large momentum transfer Q 2 between the leptons but 
a small momentum transfer (t) between the hadrons, has 
been recently demonstrated to open the possibility of ob- 
taining a quite complete picture of the hadronic structure. 
The information which can be accessed through these ex- 
periments is encoded by the Generalized Parton Distribu- 
tions, GPDs[T|[2] (for recent reviews see[3]), which give in 
particular information on the transverse location of quarks 
in the hadrons [J- Recent measurements of the azimuthal 
dependence of the beam spin asymmetry in DVCS [5 , 6 
have provided experimental evidence to support the valid- 
ity of the formalism of GPDs and the underlying QCD fac- 
torization of short-distance and long-distance dominated 
subprocesses. 

The theoretical arguments used in deriving factoriza- 
tion theorems in QCD for the nucleon[7] target case can 
be applied to the deuteron case as well, and therefore one 
can develop the formalism of GPDs for the deuteron[5]. 
From the theoretical viewpoint, it is the simplest and best 



known nuclear system and represents the most appropri- 
ate starting point to investigate hard exclusive processes 
off nuclei [S]. On the other hand, these processes could of- 
fer a new source of information about the partonic degrees 
of freedom in nuclei, complementary to the existing one 
from deep inelastic scattering. Experimentally, deuteron 
targets are quite common and as a matter of fact, DVCS 
experiments are being planned or carried out at facilities 
like CEBAF at JLab and HERMES at HERA, where some 
data have already been released [10] . One should of course 
distinguish the case where the deuteron serves merely as 
a source of slightly bound protons and neutrons from the 
case where the deuteron acts as a single hadron. In the 
former case, the scattering is incoherent and the deuteron 
will break up during the reaction. In the latter case, to 
which we devote our study, the deuteron stays intact after 
the scattering. The fact that this occurs in a non negli- 
gible fraction of events is not evident to everybody, since 
it is usual, but uncorrect, to mix the concepts of hard 
and destructive reactions. Indeed, as estimates given be- 
low will show, the very nature of deep electroproduction 
in the forward region is that the target is not violently 
shattered by the hard probe. The fragile nature of the 
deuteron thus does not prevent it from staying intact. This 
picture should of course be experimentally tested through 
the comparison of rates for coherent and incoherent elec- 
troproduction. The need for a deuteron recoil detector is 
primordial in this respect. 

The paper is organized as follows: In section 2 we 
remind the reader of the formalism of Generalized Par- 
ton Distributions for spin-1 targets in general and the 
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deuteron in particular. In section 3 we explain in detail 
the construction of the impulse approximation to evalu- 
ate the helicity amplitudes and in section 4 we derive the 
deuteron GPDs from the helicity amplitudes and study 
the properties and implications with a numerical model. 
In section 5 we give the useful formulae for calculating the 
DVCS cross section. In section 6 we show our numerical 
estimates for the usual observables in the DVCS case and 
comment on the feasibility of experiments. In section 7 
we examine the electroproduction of mesons. Througout 
the paper we will limit ourselves to the quark sector of 
the GPDs, which is a good approximation provided the 
Bjorken variable XBj is not too small. This variable is de- 
fined as usual as XBj = 2P^' i- e -> ^ is given in lab frame 

by XBj = 2A Iu , where M is the deuteron mass and v the 
virtual photon energy. Gluon effects will be needed for un- 
derstanding higher energies experiments. Previous short 
reports on our results have been presented at recent con- 
ferences [TT] 



2 GPDs in the deuteron: definitions and basic 
properties 

A parametrization of the non-perturbative matrix elements 
which determine the amplitudes in DVCS and DEMP on 
a spin-one target were given in terms of nine GPDs for 
the quark sector [5]: 
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where \P, A) represents a deuteron state of momentum P 
and polarization A, P = (P + P')/2, and n 1 * is a light- 
like vector with n + = 0, n± = 0. Due to the spin-one 
character of the target, there are more GPD's than in 
the nucleon case, but at the same time the set of po- 
larization observables which in principle could be mea- 
sured is also richer. Not much is known about these non- 
perturbative objects which encode the way quarks are con- 
fined in deuterons, except a limited set of sum rules and 
some limiting case values. Sum rules[5] relate these GPDs 
to the usual deuteron form factors : 



J dxHi(x,€,t) = Gi(t) 



dxHi(x,Z,t) = Gt(t) 



(*= 1,2,3), 



(i=l,2), (5) 



or lead to a null average : 

/ dxH 4 (x,£,t) = [ dxH 3 (x,£,t) = 0, 



D(P) 




D(P') 



N(p 2 ) 



Fig. 1. Estimating the heliticy amplitudes for the -y*D — > -yD 
in the impulse approximation. The final result is a convolution 
model between the deuteron wave function and the GPDs for 
the nucleon (upper blob). 



dxH 5 (x,£,t) = / dxH 4 (x,£,t) = 0. (6) 



Taking the forward limit of the matrix elements defin- 
ing GPDs leads to the relations [8] between GPDs and 
parton densities in the deuteron (with obvious notations) 
as : 



Hi (a;, 0,0) 



q 1 (x) + q~ 1 {x) +q°(x) 



H 5 {x,0,0) = q v {x) . 

#i(z,0,0) =ql(x)-q^{x) 



(7) 



for x > 0. The corresponding relations for x < involve 
the antiquark distributions at —x, with an overall minus 
sign in the expressions for Hi and H§. 

In all this paper we will restrict to the quark contri- 
bution. Gluon contributions are expected to be small at 
medium energies but should be included in a more com- 
plete description of the process. Therefore we will limit 
ourselves to values of ^Bj not smaller than 0.1. 



3 Helicity Amplitudes in the Impulse 
Approximation 

The impulse approximation is the zeroth order to explain 
the photon-nucleus interaction, but it is the first model 
one has to analyze since the bulk of the physics is already 
contained in it. 

In the previous section we mentioned that the rele- 
vant quantities are the deuteron GPDs . For the sake of 
simplicity we will model the matrix elements V\i\ Qand 
Aa'a CP and recover the GPDs just by using the relations 
given in the appendix. Furthermore, since we are going to 
limit ourselves to the quark content of the deuteron and 
it is an isoscalar target, we will denote: 
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and a similar relation holds for A 
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Let us denote by P p (P' M ) the momentum of the in- 
coming (outgoing) deuteron and A (A') its polarization 
state, that sometimes we will denote by 0, + or — . To per- 
form this analysis we will choose a symmetric frame where 
the average momentum P p = (P p + P ,fJ, )/2 has no trans- 
verse components. We will also need a light-like vector 
to define, together with P M , the light-cone and satisfying 
P ■ n = 1. To be more concrete, it is convenient to choose 
a frame where P M moves fast to the right. 

The momentum transfer Z\ M = P' M — P p has a longitu- 
dinal and a transverse component. The skewness controls 
the fraction of momentum transfered in the '+' direction: 



Therefore, we have that the relevant kinematical quan- 
tities for the nucleons that make up the initial deuteron 
are (a = a{) : 



Pi. 



Pi 
i4 



P2. 



a(l +i)P + 
(l-a)(l + 0P + 

2 



(17) 



and for the final deuteron we have(cv' = a^) : 



A ■ n 



(P + P')-n 2P^ 



(9) 



With these considerations, the four-vectors correspond- 
ing to each deuteron areQ: 



PM = ((1+0 p + M 2 _+A ± 2 /2 _Ax 



P>» = ((1-Z)P + , 



Af 2 + A ± 2 /2 A 



.-^) (11) 



2P+(l-£) 2 
The invariant momentum transfer is written as : 



4£ 2 M 2 



t = A A 



(12) 



The positivity of implies that there is a minimal 
momentum transfer to for a fixed £: 



to = - 



4£ 2 Af 2 



(13) 



and at the same time, for a given t, there is un upper 
bound on the allowed values for £ : 



e< 



-t 



AM 2 -t 



(14) 



The GPDs depend in addition of one more variable 
x which is defined as the fraction of average momentum 
carried by the partons in the '+' direction : 



P'i 



P? 



P 2 



= c/(l-£)P+ 

= (1 - a')(l OP" 

= 

2 



(18) 



We can now use the decomposition of the deuteron 
states in terms of nucleon states and the wave function 
defined in Appendix Al, Eq. (|7T|) . to get : 



V? 

A' A 



J da da 1 dpi ± dp' 1± J — 



i + € i 
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• @(a(l + - |a;| - 0(a(l + - 2£) (19) 

• XA'( a '> k J-> A 'l> A 2)XA(a,k^, Ai, A 2 ) 

: Ai .A2 

■ 2 / ^re irea: <^i, Ai|^ 9 (--n)T • n^ ff (-n) jpi, Ai) 

In the equation above, the variables of the spectator 
nucleon have been eliminated just by using the normaliza- 
tion properties of the one-particle states. The arguments 
of the deuteron wave function (see appendix Al for de- 
tails), k'^ and k^, are the transverse momentum of the 
active nucleon in a frame where P^_ = and Pj^ = 
respectively. Their relationship with the transverse mo- 
mentum in the symmetric frame is : 



k ■ n 
P ■ n 



(15) 



with fc M = (k^ + fc' A ')/2. Therefore, the longitudinal mo- 
mentum of the initial parton is (x + £)P + , whereas the 
final one has (x — £)P + , delivering a longitudinal transfer 
A+ = -2£P+ to the deuteron. 

Now let us turn our attention to the kinematics at the 
nucleon level and let us define the fraction of longitudinal 
momentum carried by each nucleon in the deuteron as : 



Pi 
P+ 



pr+ 



(16) 



we use the following notation 
(« + ,ti",vi), with v ± = -^(v° ± ^ 3 ) 



for the four-vectors 



k ± = ki^ = p ±± + a 
k' ± = k' = p' 1± - a' 



2 

A_l 
2 



(20) 
(21) 



The Heaviside functions in Eq. (|2U|) ensure the posi- 
tivity of the 'plus' momentum carried by the nucleons and 
put a lower bound on the integration over a. The first one 
stands for processes where \x\ > £ whereas the second one 
acts when we are on the ERBL region where \x\ < £. 

Therefore, we end up with matrix elements of a non- 
local quark operator between one-nucleon states, which is 
parameterized in terms of nucleon GPDs. To take advan- 
tage of usual parameterizations found in the literature, it 
is convenient to keep on working on a symmetric frame. 
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However, since nucleons carry some transverse momentum 
within the deuteron, the symmetric frame for the deuteron 
is not the symmetric frame for the active nucleon. 

By performing a transverse boost it is possible to eval- 
uate the matrix element in (|20l) in a symmetric frame, 
which has the property: 

Pi x + Pix = (22) 

where we have marked with a ~ the quantities in this 
boosted frame. Since it is a transverse boost, it does not 
change the '+' components of the vectors. The light-like 
vector is not changed either, since it is a light-like vec- 
tor with no transverse components, i.e. = n 11 . 

In that frame, the parameterization of the nucleon 
GPDs is made with variables that refer to the initial and 
final nucleon, i.e., we define : 



k ■ 



x N 



IN 



p 1 -h a(l + £)-£ 
A ■ n £ 



2p 1 ■ n a(l + - € 



(23) 
(24) 



Due to the lower bound on the values of a, we have 
£m > £. Moreover, it can be checked that § = f^, which 
is consistent with the fact that we are probing qq distribu- 
tion amplitudes in the deuteron only through the nucleon. 
In other words, when we enter the ERBL region in the 
deuteron (i.e. when \x\ — £), we do so at the nucleon level 
(i.e. when \xn\ — £n)- 

The transverse momentum of the nucleon that inter- 
acts with the photon is, after the boost : 



Pix 

PL 
A i 



2 

Ai 

2 

(l + 6v)A ± 
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(25) 

(26) 
(27) 



Now we can use the parameterization of the nucleon 
matrix element given in the appendix, and after some al- 
gebra and changes in the integration variables we reach 
the final result for V\<\ : 



V q - 
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and a similar expression for A q x , x 
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2Ai£jv ^° t V\iE 1S (x N ,£ N ,t)dx> i _ Al 



2m 



(29) 



The factor 2 in front of the formulae above stands 
for the number of nucleons, so that the isoscalar nucleon 
GPDs (H ls , E IS : . . . ) is the isoscalar combination within 
one single nucleon : 



H IS (x N ,£ N ,t) = i [H u (x N ,£ N ,t) + H d (x N ,£ N ,t)] . 

(30) 

The phase that goes with the nucleon helicity flip GPDs 
is given by : 



2XA X - iA y 

iXTi 



(31) 



For the sake of clarity we have omitted the Heaviside 
function in the integrals above but recall that there is a 
lower bound on the value of a, which is 



2£ \x\+£ 



i + r i+€ 



(32) 



4 Deuterons GPDs 



Once we have obtained the helicity amplitudes it is straight- 
forward to get from them the deuteron GPDs, just from 
the definitions given in Eq. (|3l4p . To do so in a simple way, 
one can use the light-cone polarization vectors given in 
[5]. The analytical expressions are summarized in the ap- 
pendix. 

At this point one may argue that definitions of deuteron 
GPDs were not actually neccessary to reach the results of 
the preceding section and that one can derive the cross 
sections and observables, directly from Eqs. (|2"5)) and (|2"5|) . 
But it should be emphasized that the genuine objects that 
parametrize the hadronic structure are the GPDs, the rest 
being just kinematics. The GPDs have well defined proper- 
ties in some limits and their analysis could help us in test- 
ing the soundness of a model, as a complementary check 
to the comparison with experimental data 
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Fig. 2. Deuteron generalized parton distributions in the im- 
pulse (a and b) and beyond (c) the impulse approximation 

4.1 Deuteron wave function 

We need a specific model for the spatial deuteron wave 
function. As can be seen in the appendix, in the lower 
Fock-space approximation one can link the light-cone wave 
function to the usual instant-form wave function through a 
identification of variables. We have chosen a parametriza- 
tion of the spatial wave function given by the Paris Poten- 
tial |12) which has a S-wave supplemented with a D-wave 
component with a probability of 5.8 %. We do not ex- 
pect a strong dependence on the chosen parametrization 
for the deuteron wave function. Most of them are identi- 
cal in the low-momentum region since they are strongly 
constrained by the well known form factors. Differences 
between parametrizations are significant only in the large 
momentum region. Since we are going to limit ourselves 
to the low-momentum transfer region, we will not be es- 
pecially sensitive to the tail of the wave function. 

Nonetheless, before going through the details of the 
results, let us discuss some features of the deuteron GPDs 
that may be expected from quite general grounds. The 
skewness parameter £ determines the momentum transfer 
in the longitudinal direction: 

A+ = (P'+ - P+) = -2£P+ , (33) 

and in the generalized Bjorken limit this is entirely fixed 
by the kinematics of the virtual photon (£ xbj/2). In 
the impulse approximation, this momentum transfer has 
to be provided by the active nucleon, and after that, the 
final state of this active nucleon still has to fit into the final 
deuteron. Since the deuteron is a loosely bound system, 
one cannot have a very asymmetrical sharing of longitu- 
dinal momentum between the nucleons and one may thus 
guess that the formation of the coherent final state will 
be strongly suppressed in the impulse approximation for 
large skewness. 

To be more quantitative let us define the longitudinal 
momentum distribution of the nucleon in the deuteron as: 

n\(a)= J \x\(/3,k A _,X 1 ,X2)\ 2 5(a~ /3) , 

(34) 

which is normalized according to 

J da n\{a) = 1 . (35) 

In Fig. [3]a we show no (a) evaluated with the wave 
function from the Paris potential [T^] . This distribution is 
strongly peaked at a = 0.5 and its width is of the order 




Fig. 3. (a) Longitudinal momentum distribution of the nu- 
cleon within the deuteron. (b) Gap between the fractions of 
longitudinal momentum carried by the active nucleon before 
and after the interaction as a function of ^(^Bj) an d a 



of the ratio of the binding energy divided by the nucleon 
mass. 

In the impulse approximation, the active nucleon af- 
ter the interaction with the photon carries a fraction of 
longitudinal momentum which is given by 

a ' = a °^—(l-a) . (36) 

In Fig. [21b we plot the difference a — a' as a function 
of a and for several values of the skewness. We see that 
for XBj > 0.1 this difference is larger than the width of the 
momentum distribution, and therefore, we will inevitably 
have a too fast or too slow nucleon (in the longitudinal 
direction). In this case the central region of momentum, 
where a maximal contribution is expected, is missed and 
then the cross sections will decrease very fast with XBj • In 
other words, there is an increasing difficulty in forming a 
coherent final state as the longitudinal momentum trans- 
fer, i.e. ccbj increases. In that case other coherent mecha- 
nisms, which could involve higher Fock-space components, 
will presumably become dominant. Not much is known 
about these states, but it should be emphasized that the 
suppression of the diagram of Fig. 1 occurs at XBj as low 
as 0.2, so that there is room to check the importance of 
the contribution of these 'exotic' states. 

The choice of the nucleon GPDs deserves a more de- 
tailed discussion. 



4.2 Modelling nucleon GPDs 

Let us first consider the helicity conserving nucleon GPDs. 
Following [13] we have taken a factorized ansatz for the t 
dependence of the nucleon GPDs: 

H u (x N ,^ N ,t) = h u (x N ,Z N )^F?(t) (37) 

H d (x N , 6v, t) = h d {x N , £ N )F?(t) (38) 
H\x Nl 6v, t) = h"(x N , Z N )F«(t) (39) 

and neglected the strange quark contributions H s . The 
flavour decomposition of the proton and neutron Dirac 
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form factor, for which we have taken the usual dipole pa- 
rameterizations [14], gives: 



F?(t) = 2F? 



pn 
■Ff 



For the axial form factor we have taken F q (t) = (1 — 
t2 ^- 2 with M A = 1.06 GeV [15 ]. For h q and h q we 
follow the ansatz based on double distributions: 



t/M\r 



h q (x N ,£ N ) 



1 rl-x' 

dx' / dy' 

ii J-l+x' 



5(x N - x' -^Ny')q(x') 



S(x N + x' -£,Ny')q{x') 



7r(x',y'), (42) 



h q (x N ,£ N ) 
w(x',y') 



/X—x' 
dy' S(x N - x' - £ N y') 
-l+x' 

Aq v (x')Tr(x',y'), (43) 



3 [l-x'f 

4 (l-x') 3 



(44) 



where we have only considered the polarization of the va- 
lence quarks. To avoid numerical problems with the inte- 
grals in the low x region we have followed the procedure 
explained in [5] . Throughout this work we have taken the 
parameterization MRST 2001 NLO [H] for the unpolar- 
ized parton distributions and the parameterization LSS01 
[T?] for the polarized ones. 

Concerning the helicity flip nucleon GPDs, E q and 
E q , we can safely neglect the latter since we deal with 
an isoscalar target. The former is suppressed in the Vy\ 
amplitudes, Eq. ([25)1 by kinematical factors. However, one 
might think that there could be physical situations which 
could be sensitive to this GPD: in the amplitudes where 
A' 7^ A the transition with H IS is done at the cost of using 
the D-wave of the deuteron, i.e., by making use of angular 
momentum. The term with E IS could flip the helicity at 
the nucleon level and therefore, the (rather small) D-wave 
admixture in the deuteron is not necessary. 

Unfortunately, most of the observables are dominated 
by amplitudes with A' = A. We have checked that effects 
due to E 1 are negligible and we have only shown how tiny 
they are for the sum rules, just for illustrative purposes. 
Furthermore, when modelling E q following the steps ex- 
plained in [13) , one realizes that the isoscalar combination 
is suppressed. Recall that E q is normalized to the Pauli 
form factor, that in the forward limit gives just the anoma- 
lous magnetic moment, very small for the isoscalar case. 

Let us stress that the available models of GPDs are 
fraught with uncertainties, in particular in the ERBL re- 
gion. There, GPDs describe the emission of a qq pair from 
the target, and an ansatz only using the information from 
usual parton densities should be used with care. Notice 
also that, while for x > £ GPDs are bounded from above 
[15] . no analogous constraints are known in the ERBL 



region. A particular type of contribution in the ERBL re- 
gion is the Polyakov- Weiss D-term [19] , which we will not 
include in our analysis. 



(40) 

M-n 4.3 Results 



With the ingredients mentioned before we plot in figures [4] 
and [5] the corresponding generalized quark distributions, 
which is flavorblind for the deuteron case. The support 
of these functions is— 1 < x < — 1 but we have plotted 
only the central region. In addition, due to the assump- 
tion made when modelling H for the nucleon (the non- 
contribution of the polarized sea) we have that Hi(x < 
£,£,i) vanishes. 

The rapid falloff of the GPDs with x reflect the fact 
that the impulse approximation, i.e, the single nucleon 
contribution cannot account for very large longitudinal 
momentum. 

Notice also the huge differences in the scales of the var- 
ious GPDs : for the vector sector, H3 dominates over the 
others, whereas H4 or H§ are very small. The respective 
sizes may be related to the values of the different deuteron 
form factors for the GPDs that have a sum-rule connection 
to them (see Eq[5]). 

The form factors that we have used in the current 
(Gi, G2, G3) are related with the usual charge monopole, 
Gc, magnetic dipole, Gm and charge quadrupole, Gq in 
the following way : 



Gi(i) = G c (t) - -r]G Q {t) 
G 2 (t) = G M (t) (45) 
(1 + v )G 3 (t) = G M (t) - G c (t) + + G Q (t) 

whith i] = jtjt. With the flavour decomposition of the 
form factors for the deuteron, we have : Gf = Gf = G\ = 
3Gi. The dominant form factor is G\ due mainly to the 
size of Gq(€) (see [20), and if we consider the form factors 
as a normalization condition for the GPDs, it is natural 
that H3 dominates over the other GPDs. 

Notice, however, that the fact that a GPD is large 
does not mean necessarily that it plays a major role in the 
observables : it has to be multiplied by the corresponding 
kinematical coefficients. 

It is worthwhile mentioning that we have plotted the 
GPDs at a particular value of —t, i.e. we cannot set t = 
to study the £ and x behaviour. The reason is that, even 
if we assume a factorized form for the t dependence in the 
nucleon, in the deuteron we cannot isolate this t depen- 
dence. In fact, there are two sources of t dependence in 
Hi and Hi : first the explicit t dependence in the nucleon 
GPDs and, the most important one, the transverse mo- 
mentum in the deuteron wave function. Then, we cannot 
circumvent the kinematical relationship between £ and t, 
Eq. (|12l) and for a non-vanishing £ we have inevitably a 
non- vanishing t. 
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4. Generalized Quark Distributions for the deuteron at 
2 GeV 2 , i = 0.1 and t = -0.25 GeV 2 . 



0.2 




-0.4-0.2 0.2 0.4 



-0.4-0.2 0.2 0.4 




-0.4-0.2 0.2 

X 



0.4-0.2 0.2 0.4 

X 



Fig. 5. Generalized polarized Quark Distributions for the 
deuteron at Q 2 = 2 GeV 2 , £ = 0.1 and t = -0.25 GeV 2 . 
Hi (upper-left), H2 (upper-right), H3 (lower-left), H4, (lower- 
right), 



4.4 Sum Rules : tests and discussions 

In the impulse approximation we have retained only the 
lowest Fock-space state of the deuteron (Fig. [2J a and b). 
As we see, the qq components which are tested in the 
region |x| < £, arc considered only within the nucleon 
itself (Fig. Hb). 

We have also neglected the NN components in the 
deuteron wave function, which could give rise to diagrams 
like the one in Fig. [2Jc. When one evaluates the elastic 
deuteron form factor [21,22] this is an exact approxima- 
tion since one can always choose a frame where the mo- 
mentum transfer is purely transverse and in that case, 
A + = and no pairs can be created from or annihilated 
into a photon. 

In the DVCS and DEMP cases there is a non- vanishing 
momentum transfer in the longitudinal direction, controled 
by the skewness parameter £. Therefore, there are neces- 
sarily diagrams where the final photon goes out from the 
anhilation of, for example, a NN pair or a qq(see figure 
[2Jc). One has to include these higher Fock-space compo- 
nents to recover Lorentz invariance. 

Lorentz invariance is the physical reason why the sum 
rules ©, obtained by integrating the GPDs over x, be- 
come ^-independent. When we perform that integration 
we have matrix elements of local operators (form factors) 
that cannot depend on the artifacts of the kinematics, i.e. 
of the separation between transverse and longitudinal mo- 
mentum transfer (see [23 for a more detailed discussion 
on this point). 

We can make use of this relationship between the £- 
independence of the sum rules and the contribution of 
higher Fock-space states in the deuteron to check how ac- 
curate the impulse approximation is. We have plotted in 
figure [5] the quantities: 



J' dxH*(x,t,t) 



(46) 



For a fixed t, the functions would be constant if 
the impulse approximation was exact : straight lines in 
the figure show the 'theoretical' values of the sum rule ac- 
cording to the experimental parameterization of the form 
factors. Any residual £ dependence of is a measure 
of the importance of the higher Fock-space states that 
we have not included in our description. Looking at this 
figure and to the corresponding one for the axial case, we 
can see that this dependence is fairly mild, which indicates 
that, in the kinematical regime that we are interested in, 
the deuteron is essentially a two-nucleon statc0. As £ in- 
creases the impulse approximation would become a too 
rough approximation with respect to Lorentz invariance. 

One should distinguish between the variation in £ of 
the quantities and the particular values they take, 

which are sensitive also to the details of the employed 



2 We checked that the inclusion of the nucleon GPD E IS 
does not introduce any improvement at all: its contribution 
vanishes exactly at £ = 0, and is always small at other values 

of e 



8 



F. Cano, B. Pire: Deep Electroproduction of Photons and Mesons on the Deuteron 



model. In figure[B]we see that the points obtained with our 
calculations are quite close to the experimental parame- 
terization. Obviously, the models works better at smaller 
£, for the reasons exposed above. In fact, in figure [7J we 
show the t-dependence of the sum rules at £ = 0, where it 
is clearly seen that results agrees quite well with the exper- 
imental parameterization, when available, or with the val- 
ues imposed by time reversal or Lorentz invariance. This 
comes as no surprise since it is well known that the light- 
cone deuteron wave function is able to give the deuteron 
form factors at the momentum transfer we are working 
(see [23] for a recent review). In the context of our discus- 
sion, at £ = the pair creation or annihilation with the 
photon vanishes in the light-cone formalism. 

One final remark concerns also the subtleties of the 
light-cone formalism: nucleons are on-shell, i.e., they ver- 
ify, with the notation employed in the previous section, 
that: 



Pi 



2p1 



(47) 



But, they are off light-cone energy shell, and as a con- 
sequence, if P + = pf +pj and = pi^ + p 2i (as it is 
the case), one has that P~ Pi +P2 ■ Therefore, one has 
that the momentum transfer t defined from the deuteron 
variables does not coincide with the one defined from the 
variables of the active nucleon. Moreover, the upper limit 



over £ff is not 



~ N in nut, 4m 2_ t ■ If nevertheless one enforces £fj to 
have this upper limit, this leads only to a tiny shift in the 
value of a m i n in the integrals. From the practical point of 
view, these differences are too tiny to be seen in the nu- 
merical calculation, unless one goes to very large values of 
—t. This just reflects the fact that the off-shellness effects 
in the light-cone energy are of the order of the binding 
energy over the longitudinal momentum |25) : 



P- 



(Pi + P2 ) K 



V 
P+ 



(48) 



and in our case this is of the order of the binding energy 
of the deuteron over the center of mass energy, i.e, very 
small. 
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Fig. 6. Sum rules for the vector GPDs at Q 2 = 2 GeV 2 and 
t = -0.5 GeV 2 . Solid lines are the expected theoretical re- 
sult and points are the results obtained with our model. Filled 

TO 

points : the nucleon GPD E is not included ; Empty points : 
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where the three form factors have been measured in the 
low and medium momentum transfer ranges 20 . The lep- 
tonic tensor L^ v is given by 



5 DVCS Amplitudes and Cross Sections 

There are two processes that contribute to the deeply vir- 
tual Compton scattering amplitude of Eq. 1 under consid- 
eration. The first one is the Bethe-Heitler process where 
the outgoing photon is produced from the lepton line. Its 
amplitude (for either electrons or positrons) is given by 



31 



BH 



t 



(49) 



where e is the proton charge. The deuteron current is given 
by: 



j„ = - GxiW ■ e)2P, + G 2 (t) [(e'* • 2P)e IM + (e ■ 2P)£] 

P, 



L^ v = u(k ', h') 



1 



A' 



Y + f (F7) 7 ? (M) 

(51) 

The Bethe-Heitler process is thus completely known in 
terms of already measured form factors. 

The second process where the photon is emitted from 
the hadronic part is more interesting in terms of the study 
of the hadronic structure. It is called virtual Compton 
scattering since it can be decomposed in a 7*A -> jB 
process. Its amplitude Tycs is written as 



(52) 



G 3 (t)(e'*-2P)(e-2P)j£ 



11 - 
(50) 



where the upper sign is for electrons and the lower one for 
positrons and the function i? comes from the decompo- 
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Fig. 7. Sum rules for the vector GPDs at Q 2 = 2 GeV 2 
and fixed £ = as a function of t. Solid lines represent the 
expected theoretical values whereas points are the results of 
our evaluation. 



sition of the leptonic current in terms of the polarization 
vector of the virtual photon 



dx 



x — £ + irj x + £ — irj 



E«; 



dx 



(56) 



x — £ + irj x + £ — irj 



A' A 



l(k' : h)Yu{k,h) = ^ 



vT^ 



J?(M)e"(9,A) (53) 



with the convention 60123 = +1. 

For completeness, let us first write down the formula 
for the cross section of the Bethe Heitler process on the 
deuteron 



x 



v2 



Ei r > 



BH 



(47ra e 



f 2 



-[K A A{t) + K B B{i)] (57) 



Sometimes we will also use: 

Q 



where A, B are the elastic structure functions of the deuteron, 
which are well known in the low momentum region, and 
the kincmatical coefficients are: 



Q(h,\) 



-.n(h,X) 



(54) 



vT^ 

The photon-deuteron helicity amplitudes are defined 



as: 



Mwx',HX = e*Jq',X')e u (q,X)H^ 



(55) 



K-B 

JCa 



2M 2 



[(2(k-q')+t) 2 + (2(k-q')+Q 2 ) 2 } 



(k ■ q')(k' ■ q' 
-Kb + „ At ,A M 2 + s + Q 2 - 2s kp ) 



{k' ■ q' 



and 



t 



- {(2M 2 + Q 2 - 2s kp f - (Q 2 - 2tf 

(58) 



2{k ■ q')(k' ■ q>) 

+4(Q 2 + s - s kp ) 2 + 4t(t + s - s kp )} 
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with Sfcp = (k +p) 2 . 

The VCS amplitude gives a contribution to the cross 
section which may be decomposed in terms of its azimuthal 
dependence as 



El* 



VCS I 



1 (47ra cm ) £ 
3Q 2 (l-e; 



53 (2\M W i,m\ 2 + 2\M H , 1 , H ^ 



H.H' 



+ 4e\Mffn H o\ 



+ Ay/e(l + e) cos(f>Re{M H %HoM^ tl „_ x - M wl ^ M* H , 1Hi 



- 4ecos(2$ Re[M H >i,H-iM£, hm ]) . 

The interference between the two processes leads to 
a contribution to the DVCS cross section which may be 
written as 



E PvcsTbh + Ty CS T B n) = To 



2 (47ro; e 



3Qt % /T — i 
2hHe[e;(q',\' = +l)j u L^{2(h,\)} 



H,H',\,h 

Re[M H , 1MX ] 



(60) 



where the upper sign stands for electrons and the lower 
one for positrons. 

The study of the initial electron helicity dependence 
may be expressed through the following weighted contri- 
butions to the cross section 



BH | 



= 



VCS I 



4 (47ra em ) 3 



(61) 



(62) 



3 Q 2 \ 1-e 



2 (47ra e 



5> (TvcsT^h + T^ cs Tbh) - T 3 

^ 2hlm[e;(q',X' = +l)j v L^"Q{h, A)] 

H,H',\,h 

Im[M H 'i,Hx] 



(63) 



(59>" 



where the upper sign stands for electrons and the lower 
one for positrons. 



6 Numerical Results for DVCS 

Our model enables us now to estimate the cross section of 
coherent deeply virtual Compton scattering on the deuteron. 
We are particularly interested by the forthcoming exper- 
iments at JLab and Hermes at DESY, and we thus shall 
present results for the kinematics of these experimental 
set ups. 
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Fig. 9. Unpolarized Differential Cross section for DVCS for 
typical kinematics at JLab (left panel) and HERMES (right 
panel). Dashed-dotted line: BH only; dashed line: DVCS only; 
full line: BH + DVCS + Interference. 



We present in Fig. [9] the unpolarized cross section for 
low (left panel) and medium (right panel) energy reac- 
tions. The Bethe-Heitler and VCS contributions are shown 
as well as their interference. The relative importance of 
these contributions depend much on the production angle 
of the final photon, as can be read from the figure. To 
discuss the feasibility of the experiment, a comparison to 
the proton target case is welcome. This is shown on Fig 
[10] for medium energy reactions. 

Coherent deep VCS is certainly not a negligible effect 
at small values of t and we can expect that this process 
should soon become observable so that some knowledge 
of the deuteron GPDs will become accessible. More than 
testing the validity of the impulse approximation, the goal 
of such an experiment is to observe some definite devia- 
tion from the impulse approximation predictions, thereby 
indicating some non-trivial short distance content of the 
deuteron. To scrutinize such effects, it is interesting to 
turn to some more specific observables, such as spin and 
charge asymmetries. 

The beam spin asymmetry is defined as 



A LU {<P) = 



d^{<j3)~da^{4>) 
da^{<j))+do±{(j)) 



(64) 



where (j> is the angle between the lepton and hadron scat- 
tering planes. The numerator is proportional to the in- 
terference between the Bethe-Heitler and the VCS ampli- 
tudes. 

Our predictions calculated with our modelized deuteron 
GPD's are shown on Fig. [TTJfor JLab and Hermes ener- 
gies. The sign of the asymmetry is reversed for a positron 
beam. Such a sizable asymetry should be quite easily mea- 
sured. It will constitute a crucial test of the validity of any 
model. 

It has been shown[5] that, asymptotically, the beam 
spin asymmetry exhibits a sin(<^>) azimuthal dependence. 
We have performed a Fourier decomposition of the asym- 
metry Am obtained for the deuteron and, indeed, we 
have checked the dominance of the sin(0) component, even 



F. Cano, B. Pire: Deep Electroproduction of Photons and Mesons on the Deuteron 



11 




Fig. 10. Comparison of the DVCS cross sections on a proton 
(dashed line) and deuteron (solid line) target for x^j =0.1 in 
the deuteron case (resp. 0.2 for the proton), Q 2 = 4 GeV 2 and 
E e + = 27 GeV. 
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Fig. 11. Azimuthal dependence of the Beam Spin Asymmetry 
as defined in the text. Left: x Bj = 0.2, Q 2 = 2 GeV 2 and 
E e = 6 GeV. Right: x Bj = 0.1, Q 2 = 4 GeV 2 and E e + = 27 
GeV. In both cases t is fixed to —0.3 GeV 2 . 
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Fig. 12. Coefficients of the Fourier decomposition of the beam 
spin asymmetry : Alu — ao + si sin (j> + S2 sin 2<f> as a function 
of t. Values of x B] , Q 2 and E\ as in Fig. [11] 



at relatively low values of Q 2 . But we also have a size- 
able sin(20) component (see figure ©), which is less sup- 
pressed than in the nucleon case. This is likely to come 
from the following fact : the Bcthe-Heitler propagators ex- 
ibit when t ^ to an azimuthal dependence in cos(^) which 
goes with terms of the order of t/Q 2 , but also with terms 
of the order m 2 /Q 2 . In the nucleon case this sin(2(f) — 
2 sin{4>) cos ((f>)- component was already seen[5] under some 
kinematical conditions, and the larger deuteron mass en- 
hances it. In that respect, we expect that the scaling regime 
signed by the dominance of the sin(0) component [2] is 
likely to be reached later in the deuteron case than in the 
proton case. 

The beam charge asymmetry 



(65) 



is also proportionnal to the interference of the Bethe- 
Hcitler and the VCS processes. Its characteristic azimuthal 
dependence is shown on Fi sfT51 Its size is large enough for 
a feasible experimental evaluation. 



7 Deep exclusive meson electroproduction 

Deep exclusive meson electroproduction[26 may be dis- 
cussed along the same lines as the DVCS reaction. The 
same QCD factorization property exists which allows to 
separate a short distance subprocess from long-distance 
matrix elements, provided the initial virtual photon is lon- 
gitudinally polarized. The same GPDs appear in principle 
in the deuteron sector, but selection rules select some of 
the GPDs. The meson production is described within the 
well established colinear approximation[27] through the 
introduction of a distribution amplitude <P(z,Q 2 ) which 
is generally parametrized through its asymptotic (in the 
sense of its Q 2 evolution) expression. 

<pP(z,Q 2 )=6fPz(l-z) 



for the p meson, and 

^(z,Q 2 ) 



%sf2fz{\ 



(66) 



(67) 



for the 7T meson, with f p = 216 MeV and j n = 92 MeV. 
The resulting amplitudes read (A, A' denoting deuteron 
polarizations as above): 



M{, x = -ie 



167ras 1 
9 Q 



dz 



*p(z) 
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Fig. 13. Azimuthal dependence of the beam charge asymmetry 
Ac for the kinematics of HERMES shown in previous figures. 



dx( 

o x - £ + le 



for the vector meson case, and 



^,a(68) 



167ras 1 



dx( - 



)^'.a(69) 



for the pseudoscalar meson case, where the isosinglet na- 
ture of the deuteron has been used to simplify the results. 
The factor in front of the matrix elements Vx'x and 

Ax'x come from the flavour decomposition of the p° and 
the 7T° , i.e., ^(|im) - \dd)). 

An interesting feature of meson electroproduction is 
the absence of a competing subrocess such as the Bethe- 
Hcitler process in DVCS. Rates are also higher than in the 
DVCS case by a factor of as/a em . The effective strong 
coupling constant has been taken as as — 0.56 as advo- 
cated in [26] ■ Vector meson (mostly p°) production selects 
the vector GPDs H 1 while pseudoscalar meson (mostly ir°) 
production selects the axial ones H l . Meson pair produc- 
tion, described in the formalism of the generalized distri- 
bution amplitudes [25], may also be calculated in the same 
way. 

We show on Fig. Q3] the prediction of our model for 
p° and 7T° electroproduction for Q 2 = 4 GeV 2 , £Bj =0.1 
and for Q 2 = 3 GeV 2 , xbj = 0.2. As in the proton case, 
the vector meson production is enhanced with respect to 
the pseudoscalar meson production. The pseudoscalar pro- 
duction is quite small since the isosinglet nature of the 
deuteron forbids any enhanced E contribution due to it 
exchange in the t-channel. 
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Fig. 14. p° (left) and tv° (right) virtual photoproduction cross 
sections at XBj = 0.1, Q 2 = 4 GeV 2 (solid lines) and xgj = 0.2, 
Q 2 = 3 GeV 2 (dashed lines). 

8 Conclusion 

We have demonstrated that deeply virtual Compton scat- 
tering and deep exclusive meson electroproduction on the 
deuteron is a feasible and promising field of study of the 
deuteron structure. More theoretical work is obviously 
needed before one can draw definite conclusions from forth- 
coming data. In particular, some of the so-called higher 
twist terms arc nccdcd|29 and should be estimated, at 
least those coming from target mass effects [30] that one 
may expect to be non-negligible. A first estimate of this 
effect has been given in the second reference of Ref. 9. The 
inverse process, photoproduction of a lepton pair, some- 
times called timelike Compton scattering is also feasible, 
and should give complementary information as was shown 
in Ref. [31] for the nucleon case. Double deeply virtual 
Compton scattering [35] may also be studied along the 
same lines. We expect that a deepened understanding of 
the short distance structure of the deuteron will emerge 
from these studies. 

We acknowledge useful discussions with M. Diehl, M. 
Gargon and F. Sabatie. Special thanks to P.A.M. Guichon 
and M. Vanderhaeghen for their help with technical calcu- 
lations. This work was supported by the EC-IHP Network 
ESOP, Contract HPRN-CT-2000-00130. 



A Appendix 

A.l Light-cone Deuteron Wave Function 

We have chosen a covariant normalization for the one- 
particle states on the light cone: 

(p' + ,p ± ',X'\p+,p u X) = (2TT) 3 2p+6(p' + -p + ) 

• 5 (2) (P±' -Px)5xx> (70) 

where p + is defined in terms of the ordinary vector com- 
ponents as p + — 771 (p° + P 3 ) and p± corresponds to the 
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components in the transverse direction. We have to evalu- 
ate matrix elements of a quark operator between deuteron 
states but eventually we have to deal with two quasi free 
nucleon states. The definition of the state of a deuteron 
with momentum P and polarization A in terms of two- 
nucleon state is: 



k z = M (£ - 1/2) 
where Mq is the free mass operator: 



A/n 



fel 



C(i-0 



(78) 



(79) 



s(i -6-6 



d Pl± dp 2± ^ 2 >(P ± 
\pf,Pi±,M;pt,P2 1 



,A 2 ) 



P2j 



then the eigenvalue equation fulfilled by X\(£> kx\ Ai, A2) 
can be interpreted as a Schrodinger equation and there- 
fore, it can be related to the wave function obtained from 
phenomenological potentials in the instant-form formal- 
^ 2 lsm. More explicitely: 
(71) 



where & = -pp- and the transverse momenta ki ± are the 
transverse momenta of the nucleons in a frame where 
P±_ = (hadron frame). We have also kept the notation 
Pi ± to make explicit the fact that the states are defined 
in an arbitrary frame but the wave function always refer 
to the hadron frame. The relationship between both set 
of coordinates is given by 



hi 



Pi± - 6-Pj 



(72) 



Let us define the measures that take part in the inte- 
grals as: 



X\(£,,k±;X 1 ,\ 2 ) = ^2 



M 



461 - 6 



1/2 



• <A 2 |ii] M -(l-^-fex)|/i2) xl(k; Mi, Ma) 



(Ai|i4(£,fc ± )|Mi) 



(80) 



The global factor in the r.h.s. of the equation above is 
just the Jacobian of the transformation from the variables 
{£, k±} to {k}. We also have the matrix elements of the 
Melosh rotation, which relate the spin in the light-front 
with the spin in the instant-form of the dynamics. Fi- 
nally, the (canonical) deuteron wave function X\(k'> Mi> M2) 
is written as 1341: 



m = d6<*6 6(1 - 6 - 6) (73) 
[dp a] = d Pl± dp 2± <S (2) (P x - Pu - Pax) ( 74 ) 

We can also perform the integrals over the transverse 
momenta in the hadron frame with the measure: 

[dk±] = dk 1± dk 2± 5^ (k 1± + k 2± ) (75) 

The wave function of the deuteron with polarization A is 
Xa(6>^1xj Ai;6)^2xj ^a) with A^ being the polarization 
of the nucleons. It is normalized according to 




By taking advantage of the properties 6 = £ = 1 — 
6 and ki ± = k± = —k 2± we can further simplify the 
notation in the wave function: 



XA(6,feix)^i;6,feax)^a) = Xx(6kx;Ai,A 2 ) (77) 

A last remark concerns the connection between the 
light-cone wave function of the deuteron and the ordinary 
(instant-form) relativistic wave function obtained with dif- 
ferent phenomenological potentials. Whereas the first one 
is expressed in terms of light-cone coordinates, the latter 
one is a function of the ordinary three-vectors ki and ful- 
fills, in general, a Schrodinger type equation. It can be 
shown [33] that, if we define the longitudinal momentum 
from the light-cone coordinates as: 



^(fe;Mi,Ma) = (16 7 r 3 ) 1 / 2 £ (_ - 1 |^ lM2 m.) 

L,mj, ,m s 

{Lll\m L m s X)Y L>mL {k)u L {k) (81) 
A. 2 Helicity amplitudes and GPDs 

We give here the kinematical coefficients that relate the 
helicity amplitudes evaluated with the light-cone helicity 
vectors and the generalized parton distributions. To be 
consistent with the choice of a right-handed set of polar- 
ization vectors made in [8], these expressions should be 
used only if A x < in the explicity evaluation of Vx^x or 
A\',\. 

By applying parity properties we can reduce the num- 
ber of independent helicity amplitudes Vx> ,x down to five 
(that we have chosen to be V++, Vqq, V- + , Vo+, V+o) 
and to four in the pseudovector case Ay,x (we keep ^4++, 

A |_, Aq + , A+o). We recall these parity and time reversal 

properties: 

Vy\ = (-1) A '- A V^ A <-A (82) 
A yx = -(-l) A '- A A_v_ A (83) 

Vva(0 = (-1) v-a Vxv(-0 (84) 
A yx (0 = (-lf- X A X y(-0 (85) 

(86) 

We denote: 
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A', A 



(87) 



A'. A 



where the sum covers only the helicity amplitudes that we 
have chosen as independent ones and the non- vanishing 
coefficients are: 



„oo _ 
c 1 — 



1 



3(1 -e 2 ) 2 
i 



(3£ 4 - 7£ 2 - 2D{1 - £ 2 ) + 2) 



1 3(1 



c l — 



3L>( 



+ - 10f + - 2)) 



C ^=3(l^Vm) P(1 " a+0 

4° - -c? + (e -»• -0 (89) 

2 



1-e 2 



D(l-a 



1 l + £ 



r o+ = _ 

2 " 2I>(l-0 

4° = -^(^ -0 

Cs - D 



r ++ - 



2£ 



r~ + - - 



2C 



i i+e 



r 0+ - 

4 1-^V 2^(1-0 
c 4 +0 = c° + (^-0 

1 



r++ - - 



„oo _ 



(l-£ 2 ) 2 

1 



^ 5 (1-e 2 ) 
1 



c R = - 



£>(1 - C 2 ) 3 



(£ 2 + 2D(1 - r) + 1) 



(2C 2 + i?(l-0) 



c ° + ^(i^V^) w ^ 2)+e) 

In the pseudovector case they are : 



(90) 



1 



(i+D(i-e)) 



g(i + ^( i-a) 

2£(l + I>(l-£2)) 

gr=-ci + (^-o 
1 



(93) 



4(i +J D(i-a) 

£ 2 -£(i-£ 2 ) 2 



Co — 



4De(l-e 2 )(l + ^(l-^ 2 )) 
(1 + 



4^2i?(l-e 2 )(l + ^(l-e 2 )) 
4 ° = -0 (94) 

g 

g-+ = 1 

3 ~ 4D(i-aa+£(i-a) 
(i+o 



g0+ = 

%° = % + {Z^-Q (95) 



i-e 2 



i 



! (i - o y/2D{\ - e) 
ci° = -ci + (^ -o 



(96) 



In the forward limit (£ — >• 0, i — > 0) we get the simpli- 
fied expressions: 



(91) 



Hi (x, 0,0) - i(2y ++ + V o) 
tf 5 (z,0,0) - (-V ++ + V 00 ) 
Hi (a:, 0,0) = A++ 



(97) 

(98) 
(99) 



(92) 



The integral over H\{x, 0,0) is the parton number, 
which imposes a serious check on the contruction of the 
helicity amplitudes in the forward limit. 



A. 3 Parameterization of the nucleon matrix elements 

In a symmetric frame, where P M = (p M +p' M )/2 has no 
transverse momentum we take the following parameteri- 
zation: 
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g e -(p',A'|^(-|) 7 + ^(^)b,A) 



u(p',X') 
dn 



a +a A c 
2m 



(100) 



u(p, A) 



2tt 
u(p',\') 



e™° (p\ A'|^(--h + 75 ^( g )\P, A) = (101) 



7+75^(^^0 + 



75 ^ 

2m 



u(p, A) 



where the integration path must be understood along the 
'-' direction. By using light-cone helicity states, which is 
close to the usual helicity in frames where the particle 
moves fast to the right and allows to get compact and el- 
egant expressions, we have to insert in the previous equa- 
tions the following results: 



«(P',A')7 + «(P, A) = 2p + ^l- ehx> 
u(p',X'h + l5 u(p,X) = 2\2p+^/l-eSxx> 

u(p',\')t^^u(p,\)=2p+ 

( e 



(102) 



A 



ow H ^ V\°x,-x 

2m 



u(p',\')— l5 u(p,\) = 2\2pi 



2m 



e 



where 



2\A X 



2m 



iA„ 



|Ax| 



(103) 



and A^ =p'^-p^. 
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